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The behavior of the main solar semidiurnal tidal mode in a 

dissipative atmosphere is studied both in a rotating spherical atmosphere 

and by means of the equivalent gravity mode approximation. The former 

involves the numerical solution of a two dimensional partial differential 

equation which (due to the presence of friction) is non- separable. The 

latter involves approximating the tidal mode at the equator by means of 

an internal gravity wave on a non- rotating plane; this approximation has 

been used extensively in earlier studies of the behavior of atmospheric 

tides in the thermosphere where viscosity assumes dominant importance. 

In the present study, dissipation is modelled by Newtonian cooling and 

Rayleigh friction, both of which are taken to increase inversely with 

mean density. Coefficients are chosen to crudely simulate the effects of 

/ • 

molecular viscosity and conducitivity. The results of this study provide' 
an opportunity to evaluate the equivalent gravity mode formalism. Our 
main findings are: 

¥ 

i) Below 130 km, where friction is unimportant, equivalent gravity 
mode results are, for all practical purposes, identical to those at the 
equator obtained from a spherical calculation. 

ii) Above 130 km amplitudes over the equator obtained from the 
spherical calculation are about 30% smaller than those obtained from 

the equivalent gravity mode calculations. Also, there is a 15° (1/2 hour) 

\ 

difference in phase. 



iii) The amplitude reduction over the equator, cited above, is 
associated with a broadening of the latitude distribution of amplitude for 
the oscillatory pressure and temperature fields within the thermosphere. 
There is also a significant variation of phase with latitude within the 
thermosphere. Associated with the above variations are significant 
changes in the latitude distribution of horizontal velocity within the 
thermosphere. 



Equivalent Gravity Modes -- an interim evaluation 

* 

1. Introduction. 

The solution of the problem of linearized internal waves (including 
tides) in an otherwise static atmosphere with friction on a rotating spher 

I 

is rendered especially difficult, because the joint presence of friction, 
rotation and sphericity lead to a mathematical problem in which altitude 
and latitude dependence are no longer separable (Chapfnan and Lindzen, 
1970). Furthermore, when the friction is due .in significant measure to 
the presence of molecular viscosity (as is the case in the thermosphere) 
the resulting equation is of eighth order in altitude. Yet, the solution of 
this problem is of substantial importance since within the thermosphere, 
where tides virtually dominate the meteorology, we have a transition 

from essentially inviscid solutions to solutions dominated by friction; 
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the behavior in the latter region is, significantly, determined by the 
complicated transition region. 

In order to obtain some insight into the above problem without 
tackling the non-separability I developed what I called the equivalent 
gravity mode formalism (Lindzen, 1970; 1971; Lindzen and Blake, 1971). 
This formalism exploited the following facts: 

i) In the absence of friction the linearized equations for internal 
waves in an otherwise static atmosphere on a rotating sphere are 
separable in their horizontal and vertical dependences; the same is 
true on a rotating (or non- rotating) plane. The horizontal geometry 



and the rotation affect the separation constants (known as equivalent 
depths, viz. , Lindzcn, 197p) which in turn determine the vertical 
structure. However, fox - the same equivalent depths one will get the 
same vertical structure regardless of the horizontal geometry or 
rotation. 

ii) As a corollary to the above consider a planar, nonrotating 
frictionless atmosphere where quite arbitrarily we identify one 
horizontal direction with the North-South direction and the other with 
the East-West direction. For any internal mode (whether tidal or not) 
in a spherical rotating fi'ictionless atmosphere, there exists a mode in 
the nonrotating planar atmosphere which, through appropriate choice 

of North-South wavenumber, will have the same East- West wavenumber 
obtaining on the sphere at the equator, as well as the same period and 
vertical structure as the spherical mode. 

iii) The effects of friction are to a large extent determined by the ratio 
of dissipative time scale to wave period. 

iv) On a. nonrotating plane the frictional problem is separable. 

The equivalent gravity mode formalism consists in taking a wave 
as described in item (ii) above and matching it to a tidal mode in a 
rotating spherical atmosphere. One then examines the behavior of the 
planar mode in an atmosphere whose vertical distributions of mean 
temperature and friction (viscosity and ion drag) match those of the 
earth's atmosphere. It was argued that the results should approximate 
thosp obtained at the equator'in a full calculation on a rotating sphe-rc. 



The problem for a rotating sphere results from the fact that 
latitudinal structure of a given tidal mode will change with altitude when 
friction becomes comparably important' with the cofiolis force; the 
changing latitudinal structure can then produce modifications in 
vertical structure. It was felt, however, that in the neighborhood of 
the equator, where Coriolis terms are zero, such effects ought to be 
relatively small. Jt was shown, moreover, that oscillations within the 
thermosphere, excited by heating within the thermosphere, were not 
very sensitive to the latitude distribution of the excitation. This last 
finding was particularly important for diurnal oscillations in the 
thermosphere which it was shown would be almost entirely excited in 
situ. Indeed the theoretical results obtained’ for diurnal oscillations 
(Lindzen, 1971, PJT-hs-) are reasonably compatible with observations of 
thermospheric daily variations. 

However, the equivalent gravity mode formalism also suggests 
that there should be, within the thermosphere, a large semidiurnal 
oscillation (temperature amplitude 200 K) which propagates into the 
thermosphere from the mesosphere where it is excited by ultraviolet 
radiation absorbed by ozone. Semidiurnal oscillations of the predicted 
magnitude have not, as yet, been found in the available thermospheric 
data, and while the data are by no means certain, there is nonetheless 
good reason to reexamine the theory. One of the first points one might 
question (though hardly the only one) is the equivalent gravity mode . 



formalism, itself. Fortunately., it is at least qualitatively possible to 
do this without integrating an eighth order non-separablc partial 

*> 

differential equation. The point is that all our arguments concerning the 

equivalent gravity mode formalism should apply as well to an atmosphere 

where the friction is due to Rayleigh friction (linear drag) with a 

rate coefficient inversely proportional to density as to an atmosphere 

where friction is due to molecular viscosity. In the former case, we 

have only to deal with a second order non-separable partial differential 

equation whose numerical integration is economically feasible. It 

thus becomes possible to compare an equivalent gravity mode solution 

with an accurate numerical solution on a rotating sphere. Moreover, 

p , 

the height dependence of the quantity ~r (where op is the 

pressure oscillation and p is the mean pressure) in an atmosphere 

with Rayleigh friction (and/or Newtonian cooling) increasing as 1/ p Q 

. / 

(where p is the mean density) is similar to the behavior of the 

oscillatory •temperature and horizontal velocity fields in an atmosphere 

with molecular viscosity and thermal conductivity (Lindzen, 1968). 

Sp 

A study, therefore, of the behavior of in the simple system 

may offer quantitative insight into the behavior of the more realistic 
system. 

.The study of a system with Rayleigh friction and Newtonian 


cooling is the purpose of this paper. 
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Z. Equations. 

For the inviscid tidal problem (see Chapman and Lindzen (1970) foi 
derivations) one may reduce the linearized equations of motion to a single 
equation in a single unknown. A convenient choice for the unknown is 


Gr - - 


j p 


where 


*Y = C i’/d v = 


T)t 9't J. 2 


t = time, z = altitude, p = pressure, p <( = mean pressure an 
tidal pressure perturbation. The equation for G is 


d = 


* ~ F (4^ r hb\(3r - — 


where CT is the tidal heating per unit mass and time, 

rr_ _L_ 9. /stn© i f ^ £l±£2^ m + j£S\ 
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f - <r^ c ’ o is the eolatitude, s the zonal number of wave, 

K,~ - , g = acceleration of gravity, H = RT^/g = local scalo height, 

o Y tidal frequency, w = earth's rotation rate and r = earth’s radius, and 

T = basic temperature. Equ. (2) may be solved by separation of variables. 

o 

The resulting equations for latitude dependence (Laplace's Tidal Equation) is 
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where hn is the separation constant and %ve have assumed 

Qr — G^Cf:') G>^ 

n 

Equ. (3) with the requirement of boundedness at the poles forms an 
eigenfunction-eigenvalue problem for {©.,] and ■ 

If we expand *CT in terms of (3 ) ' s eigenfunctions (Hough functions); i. c. , 
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th ■ 


then the equation for the vertical structure of the n mode is 
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The boundary conditions for (4) are 


H f C 4 >") = 0 


at z = 0 


(derived from the requirement W =0 at z = 0), and that no energy be 
received frorn 0 Eor waves on a planar > nonrotating atmosphere^ 

(3) is replaced by a trivial equation for which one has solutions of the 
following form 
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where Vn corresponds to a zonal wavenumber, m corresponds to a 
latitude wavenumber and 




(6 


It proves useful to allow m to be non- integral. The equation for 
vertical structure remains identical to (4). This foiuns the basis for the 
equivalent gravity mode formalism. The equivalent gravity wave being 
a wave for which o is equal to that of a tidal mode, K — (the 

spherical mode's zonal wavenumber at the equator), and m is chosen 
so that h equals the equivalent depth of the corresponding spherical mode. 
The resulting gravity wave will have the same vertical structure as the 
spherical mode (Lindzen and Blake, 1971). Moreover the gravity wave fields 
arc asymptotic to the spherical solutions as one approaches the equator. 

The assumption of the equivalent gravity mode formalism is that the 

last feature remains approximately true in the presence of friction. The 
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extent to which the approximation holds, however, still has not been 
securely determined. 

To make the problem clearer we shall explicitly introduce Rayleigh 

. * 

friction and Newtonian cooling. These enter the equations of horizontal 
momentum and energy as follows: 


•|h - a *-• v ecs & = - i- fg ( jf) - v u 


(7 
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where v, a are the Rayleigh friction and Newtonian cooling coefficients, 

•» . . * 

respectively; u the northerly, v the westerly wind velocities aAd 5 T 

the temperature; is the mean density distribution. 

In terms of the geopotential (or, more precisely, ftp/p^), the 
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final equation is 
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is related to G by: 


r ._ki a. + aHv-' 51 f. . Sl.\2. A <p\+ «_a.<i!L|]> • s 

0 - a V tVt Ci* > l U ' r i<r-& )5i(;r) < 0 * H TT'-hT 


Due to the joint presence of v and to * equ. (10) is no longer sepaiable. 
It is a two-dimensional equation which must be solved numerically. 

On the other hand the cqviations for a nonrotating planar atmosphei e 
remain separable with equs. (5) and (6) still holding. Equ. (4) becomes 



Equ, (11) must also be solved numerically, but it is only one-dimensional. 

* 

Our procedure will be to focus on the main symmetric semidiurnal 
mode for which h = 7. 85 km, q — 4 *17*/ 1 day, and s = 2. v and a are 
chosen to simulate dissipative time scales due to viscosity and conductivity; 
dissipative effects arc, therefore, negligible below about 100 km. The 
explicit choice is 

* . 
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a = 


AT 

a 


r: 1.185 X 10 


-1 

see 


exp 



n? — GO 
a -b.X'i 



(12a, b, c 


' S' simulates infrared cooling (which proves to have a negligible effect 
on the main semidiurnal mode). 


Our excitation will be due to ozone and water vapor absorption, 
both of which are located below 100 km. In this region, inviscid modal 
(Hough modes) structure is meaningful, and we, in fact, consider 
excitation of the form of the main semidiurnal mode; i. e. , in an inviscid 
atmosphere only a single mode would be excited. The vertical distribution 
of excitation is given in Chapman and Lindzen (1970). An identical 
vertical distribution of excitation is used for the equivalent gravity mode 
(viz., Lindzen and Blake, 1971). For T q we use the ARDC standard 
atmosphere. The composition of the atmosphere (and hence the gas 

constants) are assumed constant. 

Both equations (10) and (11) arc solved by the method described in 

Lindzen and Kuo (1970). In both instances a resolution of . 5 km in the 


- ... - X 

vertical is used. In both cases solutions asymptote to Or ^ G 

as -> C*? • The asymptotic regime sets on by 250 km. We thercfoie 

take ~ - — at z = 300 km. For equ. (10) we find eleven points 

C r. ^ 

from pole to pole to constitute adequate resolution; 21 points altered 
eolutions by less than 5%. We will compare the spherical solution of 

at the equator with the equivalent gravity mode solution. 
Wc will also consider the deviations of the latitude structures obtained 
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on a sphere from those of inviscid Hough modes. 

. * 

» 

3. Results. 

In Figure 1 we show the distribution with height of the amplitude 
and phase of both for the e. g. m. calculation and at the 

equator for the spherical calculation. Also shown are results at the 
equator for a dissipationles.s atmosphere. The reason for showing the 
dissipationless results will be discussed later. Restricting ourselves to 
the two dissipative cases we see that there is virtually perfect agreement 
below about 80 km. Above 130 km noticeable differences begin to develop. 
Above 200 km these differences amount to about 1/2 hour (or, 
equivalently, 15°) in phase; amplitudes for calculations on a rotating sphere 
arc about 30% less than those obtained from e.g.m. calculations. The 
latter difference is not surprising: when one has a transition from 
exponentially increasing amplitudes to constant amplitude with height , ^ 


This is precisely the behavior of horizontal velocity and temperature 
when dissipation is due to molecular viscosity and conductivity. 

small variations in the height and nature of the transition can lead to 
significant amplitude differences above the transition. 

Between 80 km and 130 km there arc small differences between 
e.g.m. and spherical calculations (amounting to a few minutes of 
phase and a few percent in amplitude) which, while too small to be of 



practical significance, appear to be real. Some idea of the origin of these 
small discrepancies may be obtained by comparing e.g.m. results with 
dissipationless results in Figure la. The differences below 130 km arise 
from the fact that dissipation varying with height can cause the 
reflection of internal gravity waves. This matter is reviewed in Lindzcn 
(1970). On the other hand, below 130 km dissipationless and dissipative 
results for a rotating sphere arc virtually identical, suggesting that the 
reflectivity has been substantially diminished. 

A study of the variation of the semidiurnal tide with latitude, as 
obtained from the spherical calculation, casts some light on the 
results shown in Figure 1; it also shows some rather unexpected 
aspects of the interaction of friction and rotation in a spherical atmospher 
In Figures 2 and 3 we show, repsectively, the variation of amplitude 
and phase with height for (i* e * » the fractional variation of 

pressure due to the main semidiurnal mode) at different latitudes. 

We see in Figure 2 that the amplitude variation is similar at all’ latitudes 
for z 2 130 km -- indicative of the fact that the tide propagates 
vertically without change of horizontal form. Similarly, we see in 
Figure 3 that the phase is independent of latitude below 130 km. These 
results arc exactly what is obtained from inviscid theory when only the 
main mode is excited. However, above 130 km we see that the a-fnplitudc 
profiles draw closer together. More explicitly, I -jT“ \ at y = . 8 

(latitude = 53. 2°) increases up to 170 km, asymptoting to a constant 
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above this height; at y = 0 (equator, y - cos 0), on the other hand, ! 
increases up to 150 km, then decreases about 30% before asymptoting to a 
constant above 200 km. This leads to a latitudinal spreading of the 
semidiurnal main mode between 130 km and 200 km. As can be seen in 
Figure 3, the transition region also produces a variation of phase with 
latitude that persists throughout the upper thermosphere. This variation 
amounts to about 90° of phase between y = 0 (the equator) and y = 0. 8 
(latitude = 53. 2°). 

The alteration of latitude structure within the transition region 

(130-200 km) is the most significant result of the spherical calculation. 

Although our artificial choice of dissipative processes precludes a 

direct association of calculated magnitudes with these that might be 

expected in the thermosphere, the calculated changes of latitude structure 

should be Qualitatively indicative of what might happen in the real 

atmosphere. These changes are displayed in Figures 4a, b, c, and d 

where the variation of amplitudes (units arbitrary) and phase with latitude 

both below 100 km and above 200 km are shown for -_L , u\ v' and 

to • 

6 T respectively. In the frictionless region both ( and 

have identical horizontal structures (Hough functions). Moreover, the 
latitude distribution of both is broadened in the transition region. However, 
the variation of phase with latitude (above 200 km) is greater for 6T 

Sr 


than for 


. The effect of the transition zone on the 


latitudinal structure of (ul and (v\ is rather striking. Below 
the transition zone, the westerly velocity \ v\ , has a relative 
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minimum at the equator; 1 v\ increases away from the equator to a 
maximum near 30° which is about 17% greater than the value at the 
equator and then decreases to zero at the poles. Above the transition 
zone, ( v ( decreases monotonically away from the equator. The 
northerly velocity, which by symmetry is zero at the equator reaches a 
maximum near 40° below the transition, and hear 30° above the transition 
zone. For all fields, there are no phase variations with latitude below 
the transition zone; above the transition zone the phases of semidiurnal 
oscillations at the equator lead those at 52 by 2 to 4 hours. 


4. Implications for the thermosphere. 

The preceding results show in detail how friction and rotation 
interact in an atmosphere where friction increases exponentially from 


an inconsequential value below 130 km to values which dominate atmospheric 
behavior above 250 km. While the specific model for friction (and thermal 
dissipation) is not realistic, the results give us insight into what is 
omitted when we adopt the equivalent gravity mode formalism. The 
effect of primary importance is the broadening of the latitudinal 
distributions of fields like pressure and temperature with the con- 
commitcnt 30% decrease of oscillatory amplitudes in the equatorial 
thermosphere. Similar broadening may be expected to occur when 
friction is due to molecular viscosity. Thus in Lindzcn (1971) the 
equivalent gravity mode formalism was used to predict a thermospheric 
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semidiurnal temperature oscillation witli an amplitude of 180 K or more 
at the equator; the present calculations su 6 gcst that calculations on a 
rotating sphere might predict amplitudes on the order of only 120°K. 

As already pointed out in Lindaen (1971). the latter value is still 
large-r than what current observations suggest. 

Similarly, the implication of the present calculations that one 
shdulcl observe a phase variation with latitude should hold for more 
realistic models of dissipation as well. Such behavior has obvious 
ramifications for the analysis of data -- which in some instances 
involves averaging over latitude, a procedure, which for phase variations 
on the order of those in Figure 4d, would lead to significant underestimates 
of semidiurnal amplitudes. 

Finally, it should be added that the accuracy of the equivalent 
gravity mode formalism at the equator is within the limits suggested m 
Lindzen (1971). However, Figure 4 implies that away from the equator - 
account must be taken of the variation of latitudinal structure. 
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Figure Legends 

Fig. la. Amplitude of £il for the main solar semidiurnal 

mode at the equator for equivalent gravity mode calculation 
(with dissipation), and rotating, spherical calculations with 
and without dissipation. (N. B. 21 pts. in latitudes have been 
used for higher accuracy in the spherical calculations shown. ) 
Fig. lb. Same as la but for phase. 

Fig. 2. Amplitude of ^ for the main solar semidiurnal 

1 O 

mode at various latitudes (N. B. y = cos of latitude) for 
rotating, spherical calculation with dissipation. 

Fig. 3. Same as 2 but for phase. 

Fig. 4a. Amplitude (arbitrary units) and phase (arbitrary origin) of 

^ F as functions of latitude above 200 km (dashed) and 

Po 

below 100 km (solid). 

Fig. 4b. Same as 4a but for u'. 

Fig. 4c. Same as 4a but for v'. 


Fig. 4d, Same as 4a but for fcTV 












